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The elimination of decoherence of a multiphoton two-state quantum system by using an ap- 
propriate external driving field is considered. The multiphoton process caused by the noise field 
has a supersymmetric Lie algebraic structure. The time-evolution equations for the off-diagonal 
elements of the density operator of the two-state system are derived in the interaction picture. A 
simple explicit formula is given for the time-dependent external driving field used to eliminate the 
decoherence of the multiphoton two-state system. 
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I. INTRODUCTION 



During the past decade, quantum computation, which utilizes the superposition of many states of a quantum system, 
has attracted extensive attention [1-5]. Quantum computers can carry out computational tasks that are intractable 
for conventional computers since a quantum computing system composed of qubits (two-level quantum systems) has 
a mysterious property of quantum coherence (including the entanglement or quantum correlation), which has no 
counterpart in the classical realm [5]. However, because of the existence of the decoherence it is not clear whether we 
can finally build a truly practical and effective quantum computer [1]. Decoherence is essential to understanding how 
classical physics emerges from quantum mechanics. It is such a process by which a quantum-mechanical system that 
interacts with its natural environment loses the characteristic properties that distinguish quantum mechanics from 
classical physics. In other words, decoherence arises from the unavoidable interaction between the quantum systems 
and the noise fields, which are mostly responsible for damping of the quantum coherence of systems such as spin 
and two-state systems. Thus, if the decoherence occurs in the computation, some qubits become entangled with the 
environment and consequently the state of the quantum computer collapses [1]. Recently, much attention has been 
attracted to the theoretical and experimental work for topological quantum computation [6-12]. It is well known that 
due to its topological and global nature [13], the geometric phase of a spinning particle will not be affected by the 
random fluctuation arising in the evolution path [6,7]. This means that the geometric phase shift gate (should such 
exist) may be robust with respect to certain types of operational errors. Since in a quantum computer an error-tolerant 
quantum logic gate is particularly essential for realizing a quantum information processor, the quantum computation 
based on such a geometric phase shift gate will inevitably become an ideal scenario to this purpose [11]. Indeed, in the 
experiments performed at the end of last century, such conditional geometric phase shift gates were realized through 
the nuclear magnetic resonance (NMR) under an adiabatic condition [14-17]. The result of quantum computation 
will be exact only when the adiabatic condition can be fulfilled {i.e., only when the time-dependent parameters of 
the Hamiltonian of qubits vary extremely slowly) [7]. Once the adiabatic condition is satisfied, we meet, however, the 
problem of quantum decoherence caused by the inevitable interaction between the qubits (spin) and the environments. 
The existence of such decoherence effects demand that the quantum computation process (in which the maintenance 
of coherence over a large number of states is important) must be completed within the decoherence time. If not, the 
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error caused by the decoherence will inevitably increase. Thus, we need to change the parameters of Haniiltonian 
rapidly. However, as discussed before, such rapid changes may break down the adiabatic requirement and leads to the 
nonadiabaticity error in the results of quantum computation [7]. Therefore, it seems that there is inevitably a conflict 
between the adiabatic condition and the requirement for removing the decoherence effects. In the literature, there 
are several methods, such as error-avoiding codes [18], error-correcting codes [19] and decoherence-avoiding scheme 
[20], which can be used to reduce the decoherence effects in qubits. In the decoherence-avoiding scheme [20], the 
interaction between the two-state quantum system and the environment (such as a noise field, thermal reservoir and 
bath) is eliminated by using an external controllable driving field. In the present paper, we study the maintenance of 
coherence through the decoherence-avoiding scheme in a two-state system. Note that in order to consider the general 
interaction between the two-state system and the noise field we treat the multiphoton process, i.e., the total number 
of photons created or annihilated in each two-state transition process caused by a noise field is greater than one. 



II. INTERACTION OF A TWO-STATE SYSTEM WITH THE NOISE FIELD AND THE DRIVING FIELD 



The interaction between the two-state system and the noise field can be described by the supersymmetric multi- 
photon model [21-25], and the Hamiltonian under the rotating wave approximation can be written as [21,22] 



H 



Wo 



(1) 



where the creation and annihilation operators and a for the photons obey the commutation relation [o, a^] — 1, 
a± and (Jz denote the two-state system operators satisfying the commutation relation [cr2,cr±] = ±2cr±, and g and 
g* are the multiphoton coupling coefficients. The integer k is the rnimber of photons created and annihilated in each 
two-state transition process, and ix'o(i) andLo{i) the two-state transition frequency and the mode frequency of photons, 
respectively. Defining the following algebraic generators [23-25] 



Q = (a^fa- = 
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(2) 



one can easily show that (A'^, N ,Q,Q^) form a set of supersymmetric generators and possess the supersymmetric Lie 
algebraic structure, i.e., 



0, [N,Q]=Q, 



Q2 = (Qt)2^0, [Q\Q]=Naz 

[N,Q^]=-Ql {QIQ}=n', {Q,az}= {Q\az} =0 

[0,<T,]=2Q, [Qlaz]=-2Q\ [Q^ - qf = -N , 

where {} denotes the anticommuting bracket. 

To eliminate the quantum decoherence, an external classical field E{t) is applied to the above system. The Hamil- 
tonian for the interaction between the two-state system and the driving field E{t) is given by 



(3) 



H^ = ~{dEa++d*E*(j. 



(4) 



where d and d* denote the transition dipole moments of this two-state quantum system. The total Hamiltonian is 
-fftot = H + He- Note that here the environmental noise field (weak) that is coupled to the two-state system is 
modelled by a quantized electromagnetic field while the external driving field (strong) is expressed by a set of classical 
quantities, E and E* . 

The Hamiltonian of the above multiphoton two-state system in the interaction picture can be obtained by using 
the following unitary transformation 

i 



V{t) = exp 



(5) 



With the help of Hi{t) = V^t) (iftot - i^) V{t), one can arrive at 

Hi{t) = gexp{-i6t)Q + g* exp{iSt)Q^ - i {dEe''^''*a+ + d*E*e-^°*a-) , 
which is the Hamiltonian in the interaction picture, where the frequency detuning 5 = kuj — ujq. 



(6) 
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III. EVOLUTION OF THE OFF-DIAGONAL DENSITY MATRIX ELEMENTS 



In this section, we use the MarkofF approximation [26,27] to obtain the time-evolution equation of the off-diagonal 
elements of the density operator for this two-state quantum system drivciii by the external elassical field E{t). 

It is well known that the density operator of the two-state system undergoing a multiphoton interaction satisfies 
the following Liouville equation 

i^ = [H,{t),p,{t)]. (7) 
Let Psi(t) denote the reduced density operator of the two-state system. It follows from Eq. (7) that 

p,i{t) = TV.pi(t) = - /* Trr[Hi{t), [Hi{t'), pi{t')]]dt'. (8) 
Jo 

Here dot denotes the derivative with respect to time. If the thermal reservoir is assumed to be quite large, it does not 
change much during the time evolution process of the two-state system. Then we have pi{t) ~ psi(OPri(0), where the 
density operator of the thermal reservoir is pri(O) ~ ^^^exp{— f3Hk) /Tr exp{— (3Hk) with /3 = l/(kB.T), = Wkflt^k- 
If we take E{t) ~ E{t'), pi{t') ~ pi{t) (Markoff approximation) [26,27], Eq. (8) can be rewritten as 



Psi{t) = - f Trr[Hi{t),[Hi{t'),p,i{t)pri{0)]]dt', 
Jo 



(9) 



where the integrand can be rewritten as 



Tr,[JTi(t), [Hi(t'),Psi(i)pri(0)]] = TT,[Hiit)Hiit')psi{t)p,i{0) - i/i(f)p,i(f)prf(0)Hi(i') 

-Hi{t')psi{t)p,i{0)Hi{t) + p,i{t)p,i{0)Hi{f)Hi{t)]. (10) 

After some tedious derivations, we obtain the explicit results for the four terms on the right-hand side of Eq. (10). 
The first term can be rewritten as 

Ti,[Hi{t)Hi{t')p,i{t)pnm = < Hi{t)Hi{t')p,i{t) >r, 

= g*g {exp[7;5(t' - t)] < QQ^ p,i{t) >r + exp[-i,5(t' - t)] < Q^Qpsi{t) >r} 

+ ^d*dE*{t)E{t)exp[iiOo{t-t')]p,i{t), (11) 

where the relation = (Q^)^ = has been applied to the calculation. Here the subscript R means the average over 
the states of the thermal reservoir (noise field). The rest of the terms on the right-hand side of Eq. (10) are given by 

-Tr,[Hi{t)p,i{t)p,i{0)H,{t')] = - < H,{t)p,i{t)Hi{t') >r 
= -g^exp[-iS{t + t')] < Qp,i{t)Q >r -{g*fexp[tS{t + t')] < p,i{t)Q^ >r 
-g*g {exp[i(5(t' - t)] < Qpsi{t)Q^ >r + exp[-id{t' - t)] < Psi{t)Q >r} 

(P (d*)^ 
-—E^{t) exp[iwo{t + t')]a+p^i{t)a+ - {E*f (t) exp[-iwo{t + t')]a-p^i{t)(T- 

{E*{t)E{t) exp[iLOoit - t')]a+psi{t)a- + E*{t)E{t) exp[-itoo{t - t')]a-p,i{t)a+} , (12) 

-TT,[Hi{t')psi(t)p,i(0)Hi{t)] = -g^exp[-iS{t + t')] < Qp,i{t)Q >r -{g* f exp[i5{t + t')] < Q^p,i{t)Q^ >r 
-g*g{exp[-iS{t' - t)] < Qp,i{t)Q'' >r +exp[iS{t' - t)] < QVsi(i)Q >r} 

d^ (d*)'^ 
-—E^t) exp[iuJoit + t')]a+psi{t)<j+ - {E* f (t) exp[-iuJo{t + t')]a-Psi{t)a- 

-— {E*{t)E{t) exp[icoo{t' - t)]a+p,i{t)a- + E*{t)E{t) exp[-ja;o(i' - t)]a-p,i{t)a+} (13) 

and 
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TrMt)pri{0)Hi{t')Hi{t)] = g*g {exp[-i^(t' - t)] < Psi{t)QQ^ >r +exp[i(5(t' - t)] < Psi{t)Q^Q >k] 

+ ^E*{t)E{t) exp[-ia;o(t - t')]p,i{t). (14) 
Thus, the integrand in Eq. (9) can be rewritten in the following form 

Tr,[iri(i), [i?i(<'),Psi(i)Pri(0)]] =ri+% + % + T4, (15) 

where 

^1 = 9*9 {exp[?(5(i' - t)] < {a^fa^ >r a-a+ + cxp[-2(5(f' - t)] < a''{a^)'' >r, a+a-} psi(t) 

+g*gpsi{t) {exp[-iS{t' - t)] < (at)'^a'' >r (J_cr+ + exp[i(5(i' - t)] < a'=(a^)'= >r cr+cr_} , 
T2 = -2 exp[-i5(t + 1')] < Qpsi{t)Q >R +(<?*)' exp[i<5(t + i')] < QVsi(t)Qt >j^} ^ 
^3 = -5*5 {cxp[i5(t' - t)] + exp[-i,5(t' - t)]} (< Qpsi(t)g^ >r + < pA{t)Q >r) , 

T4 = —E*{t)E{t) {exp[iwo(t - t')] + exp[-iwo(i - t')]} Psi(i) 

-y exp[ja;o(i + t')]«^+Psi(i)'^+ - (i^*)' (*) exp[-ia;o(i + t')]a_p^i{t)c7_ 

d 

—E*{t)E{t) {exp[icjo(t - t')] + exp[-iu;o(t - t')]} [a+p,i{t)a- + (7_p,i(t)(7+] . (16) 

Assume that the reservoir is in a multiphoton state |m), where the occupation number m > k. Using the relations 
aHa^ f \m) = ^ |m), < aHa^)" >R= < {a^)V >r= < (at)'^(at)'^ >r= 0, a+|+) = = 

and = (— |cr+ = 0, one can obtain 



-I^i|+) =5*5exp[z,5(i'-i)] 
-|^2|+)=0, 



(m + A;)! ^ 



ml (m — A;)! 



-\psl{t)\+), 



(-IT3I+) = -g*g{eMiS{t' - t)] + exp[-i^(t' - t)]} i^I^±J^{+\p,,{t)\-), 
HW = ^E*{t)E{t) {eMiMt - t')] + exp[-iwo{t - t')]} (-|Psi(i)|+) 
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[E*{t)]'eM-iMt + t')](+\psim-), (17) 



where the relations {—\a+psi{t)a-\+) = {—\a-psi{t)a+\+) = has been inserted. 

Define p-+{t) = {-\p^i{t)\+) , p+-{t) = (+|psi(i)|-) and p-+{t) = {-\p,i{t)\+}. It follows from Eq. (9) that 

p_+(i) = - f TT.{-\[H,{t), [Hj{t'),p,j{t)prim]\+)dt'. (18) 
Jo 

Substituting Eqs. (15) and (17) into the above equation and neglecting the rapidly oscillating term {i.e., the last 
term) in (— 174|+), we obtain 

p.+ it) = Ci{t)p^+{t) - C2(i)p+-(t), (19) 

where the time-dependent coefficients ci {t) and C2 {t) are given by 



ci(t) = -9*9 ^ + t;;^^ i^^^ + '^,E*{t)E{t) sinu^ot, ^^^^ 



{771 — k)\ 

C2{t) = -g*g^^ ' 



1— exp( — 2<5t) . 1— exp(i<5i) 



Note that the coefficient ci{t) is complex and C2{t) has a real value. 
The complex conjugation of Eq. (19) gives 

p+_(t) = cl{t)p+^{t) - C2{t)p.+ {t). (21) 

Eqs. (19) and (21) are the time-dependent equations of the off-diagonal elements of the density operator for this 
multiphoton two-state quantum system. 
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IV. ELIMINATION OF DECOHERENCE IN THE MULTIPHOTON TWO-STATE SYSTEM 



In order to maintain the coherence of the two-state system interacting with the noise field, we should analyze Eqs. 

(19) and (21), which govern the time evolution of the off-diagonal matrix elements of the density operator, ft is 

obvious that in the presence of noise field the decoherence can be eliminated if p-\ {t) = p |_(t) — > is satisfied. This 

leads to the requirement 



CiCi 



In what follows, we will solve E{t). For convenience, we first define 



a{t) 



-9 9 



{m + k)l 



+ 



ml 



1 — exp(— ii5f) 



ml (m — 
and then rewrite the coefficient ci {t) as follows 



iS 



r,/ \ d, d . , . 

pit) = - — smtoot, x{t) 



E*{t)E{t), 



ci{t)=a{t)+l3{t)x{t). 
According to Eq. (22), one obtains {a* + f3x){a + fix) — C2 = 0, which can be rewritten as 

The two roots of the above quadratic equation are of the form 



_ g+g 



where 

a + a = —g 



x± 



(to + k)l ^ to! 



w 



\/ia-a*) 



to! (to — fc)! 



2 sin 5t 



a- a =-g g 



(TO + fc)! to! 



to! (to — fc)! 



2(1 - cosSt) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



If we apply an external time-dependent driving field with intensity given by Eq. (26), in principle the decoherence 
of the two-state quantum system interacting with the environment (noise field) can be eliminated. Fig. 1 shows the 
external time-dependent driving field required for the elimination of the decoherence of the multiphoton two-state 
system with various k (the number of photons in each transition process) and 5 — Q (without detuning) . It follows 
from expressions (20) and (27) that all the parameters C2, a + a* and a — a* are inversely proportional to the 
frequency dctiming (5. It is thus shown that according to Eqs. (19)-(21), the two-state system may not be subject to 
the quantum decoherence significantly in the case of large detuning. In what follows, as an illustrative example, we 
will consider only the maintenance of the coherence for the case of small detaining. Fig. 2 demonstrates the required 
time-dependent expression for the external driving field when the frequency detuning is (5 = O.OfcJo- It should be 
noted that in some time ranges (e.g., 2.2 x 10^^" < t < 2.5 x 10^^" s for the case of fc = f) in Fig. 2, both the 
field intensities x+{t) and X-{t) obtained in (25) are negative, i.e., E*{t)E{t) < 0, which has no physical meanings. 
So, in the time ranges where E*{t)E{t) < 0, what one should do is to set E{t) — experimentally. This, however, 
means that if we utilize the external driving field illustrated in Fig. 2, the coherence cannot be maintained completely 
after t = 2.2 x 10^^" s (for the case of fc = f). To overcome such a difficulty, here we will suggest a potential 
scheme of cyclically driving field: specifically, if we choose the evolutional behavior of the external driving field of the 
time range [0,T] (T = 2.2 x 10-^" s) and apply it to any ranges [nT, (n -\- f)r], namely, E{nT + t) = E{t), where 
t e [0,T], the quantum decoherence of the above two-state system undergoing a multiphoton interaction caused by 
the environmental noise field will be greatly reduced. 

To summarize, the multiphoton process, which possesses a supersymmetric Lie algebraic structure, is caused by the 
interaction between the two-state system and the noise field of the environment. This will lead to the decoherence of 
the two-state quantum system. In the present paper, we have considered the reduction of decoherence by using an 
appropriate external time-dependent driving field. We have used the Markoff approximation to derive the equations 
governing the time-evolution behavior of the off-diagonal elements of the two-state density operator in the interaction 
picture. A simple explicit formula for the intensity of the external driving field has been given for the elimination of 
the decoherence of the multiphoton two-state system. 

Acknowledgements This work is supported by the National Natural Science Foundation of China under Project 
Nos. 90101024 and 60378037. 
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FIGURE CAPTIONS 



Fig. 1. The external time-dependent driving field required for the elimination of the decoherence of the two-state 
system coupled to a noise field. Here m = 100, wq = 1.0 x 10^^ s~^ and the frequency detuning S vanishes. 



Fig. 2. The external timc-dcpcndcnt driving field required for the elimination of the dec;()lierence of the two-state 
system coupled to a noise field. Here m = 100, coq = 1.0 x 10^^ s~^ and the frequency detuning 5 = O.OIwq. 
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